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We develop an Eliashberg theory for multi-scale quantum criticality, considering ferromagnetic 
quantum criticality in the surface of three dimensional topological insulators. Although an analysis 
based on the random phase approximation has been performed for multi-scale quantum criticality, 
an extension to an Eliashberg framework was claimed to be far from triviality in respect that the 
self-energy correction beyond the random phase approximation, which originates from scattering 
with z = 3 longitudinal fluctuations, changes the dynamical exponent z = 2 in the transverse mode, 
explicitly demonstrated in nematic quantum criticality. A novel ingredient of the present study is to 
introduce an anomalous self-energy associated with the spin-flip channel. Such an anomalous self- 
energy turns out to be essential for self-consistency of the Eliashberg framework in the multi-scale 
quantum critical point because this off diagonal self-energy cancels the normal self-energy exactly 
in the low energy limit, preserving the dynamics of both z = 3 longitudinal and 2 = 2 transverse 
modes. This multi-scale quantum criticality is consistent with that in a perturbative analysis for the 
nematic quantum critical point, where a vertex correction in the fermion bubble diagram cancels 
a singular contribution due to the self-energy correction, maintaining the z = 2 transverse mode. 
We also claim that this off diagonal self-energy gives rise to an artificial electric field in the energy- 
momentum space in addition to the Berry curvature. We discuss the role of such an anomalous 
self-energy in the anomalous Hall conductivity. 

PACS numbers: 



I. INTRODUCTION 

Investigation on the role of the momentum-space Berry 
curvature has driven an interesting direction of research 
in modern condensed matter physics, covering weak anti- 
localization, anomalous charge and spin Hall effects, and 
topological insulators [iHll. ^ recent trend is to study 
interplay between the topological band structure and in- 
teraction, suggesting that topological excitations such as 
skyrmions and vortices carry fermion quantum numbers 
of electric charge or spin via the associated topological 
term [j]. 

Quantum criticahty has also driven an important di- 
rection of research in modern condensed matter physics, 
particulary focusing on the nature of non-Fermi hquid 
physics and mechanism of superconductivity out of the 
non- Fermi liquid state Combined with the topolog- 
ical band structure, critical degrees of freedom given by 
topological excitations in the duality picture will carry 
fermion quantum numbers, thus a novel type of non- 
Fermi liquid physics may arise. 

In this paper we study ferromagnetic quantum critical- 
ity in the surface of three dimensional topological insula- 
tors from the side of a disordered phase. A recent analy- 
sis based on the random phase approximation (RPA) has 
shown a multi-scale quantum critical point, where longi- 
tudinal spin fluctuations are described by the dynamical 
critical exponent z = 3 while transverse modes are given 
by 2; = 2 0. Although scaling can be demonstrated 



within this analysis, stability of the RPA dynamics is 
not guaranteed beyond this approximation. In order to 
describe scaling near quantum criticality, one should go 
beyond the perturbative analysis, introducing quantum 
corrections fully self-consistently. An Eliashberg theory 
is desirable, where self-energy corrections are taken into 
account fully self-consistently in the one loop level 
The scaling expression of the free energy can be derived 
within the Eliashberg approximation, where critical ex- 
ponents are determined [8'] . Unfortunately, an extension 
to an Eliashberg framework was claimed to be far from 
triviality in multi-scale quantum criticality because the 
self-energy correction, which originates from scattering 
with z = 3 critical fluctuations, turns out to change the 
dynamical exponent 2; = 2 in the transverse mode, ex- 
plicitly demonstrated in nematic quantum criticality [oj . 

In this study we construct an Eliashberg theory for 
multi-scale quantum criticality, considering ferromag- 
netic quantum criticality in the topological surface. A 
novel ingredient is to introduce an anomalous self-energy 
associated with the spin-flip channel. Such an anomalous 
self-energy turns out to be essential for self-consistency of 
the Eliashberg framework because this off diagonal self- 
energy cancels the normal self-energy exactly in the low 
energy limit, preserving the dynamics of both z — 3 and 
z — 2 critical modes. This multi-scale quantum critical- 
ity is consistent with that in a perturbative analysis of 
the nematic quantum critical point, where a vertex cor- 
rection in the fermion bubble diagram cancels a singular 
contribution due to the self-energy correction, maintain- 
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ing z = 2 for the transverse mode [9| . We claim that such 
an off-diagonal self-energy generates an artificial electric 
field in the energy-momentum space [l^ in addition to 
the Berry curvature. Based on the semi-classical approx- 
imation for the Hall conductivity, we find an anomalous 
divergent behavior at zero temperature. We discuss phys- 
ical implication of this phenomenology in connection with 
the topological nature of Dirac fermions. 

II. HERTZ-MORIYA-MILLIS THEORY IN THE 
SURFACE OF THE TOPOLOGICAL INSULATOR 

A. Model 

We start from an effective field theory for ferromag- 
netic quantum criticality in the surface of topological in- 
sulators 

£ = Vlil^r - fJ-)Saa' + iv f Z ■ (cT X V)c,a']'4>a' 

+|a,0,p + vl\v4>,\^ + mg|0,p + y(|0,|) 

+9't>i'^Waa''^a' ^ (1) 

where ip — ^ ^ represents an electron field in the 

Nambu-spinor representation and (f) = + ipyV corre- 
sponds to a spin excitation. We consider the case of an 
XY spin for simplicity, which can be generalized into the 
case of 0(3) spin, f/ is the Dirac velocity and fi is the 
chemical potential for surface electrons, cr is the Pauli 
matrix, associated with the spin state. Vb is the spin- 
fluctuation velocity and is the inverse of correlation- 
length for spin fluctuations, corresponding to an XY or- 
dered state when < and a quantum disordered 
phase when > 0. V{\4>i\) \s an effective potential 
for spin fluctuations, where an easy plane anisotropy al- 
lows the XY spin dynamics only, g is an effective cou- 
pling constant between electrons and spin excitations, 
where tuning g leads the spin correlation-length to di- 
verge, resulting in quantum criticality. The quantum 
critical point between the XY ordered and disordered 
phases is defined by the vanishing effective mass for spin 
fluctuations, tuned by the coupling parameter g. 

This effective field theory would be realized when an 
insulating ferromagnet lies on the surface of topologi- 
cal insulators. In this case ferromagnetic spin fluctua- 
tions denoted by 0^ describe those in the insulating fer- 
romagnet. Then, the coupling parameter g represents 
the strength of couplings between spins in the ferro- 
magnet and itinerant electrons on the topological sur- 
face. Another similar setting is that ferromagnetism 



will be induced by doped magnetic impurities. It was 
demonstrated that Ruderman-Kittel-Kasuya-Yosida in- 
teractions between doped magnetic impurities are ferro- 
magnetic TT] . Through this exchange mechanism, mag- 
netic atoms are expected to form a ferromagnetically or- 
dered fllm, deposited uniformly on the surface of a topo- 
logical insulator. In this case a possible ferromagnetic 
transition may be realized, controlling the distance of 
doped impurities in the regular impurity pattern. 

We see two important features in dynamics of elec- 
trons on the topological surface state. The first is that 
spin dynamics is quenched to the orbital motion of sur- 
face electrons. As a result, one finds that the role of 
an external magnetic field differs from that in graphene, 
where the pseudo-spin quantum number is constrained 
to the orbital motion. In particular, one of the authors 
could reveal that the Kondo effect and Friedel oscillation 
around the magnetic impurity should be modified in the 
presence of an external magnetic field, compared with the 
graphene structure [l2l |. The second aspect is more fun- 
damental, that is, only an odd number of Dirac fermions 
is allowed to appear in the topological surface state [2]- 
As a result, quantum anomaly is realized, giving rise to 
nontrivial topological properties of the system |13| . For 
example, a vortex in the XY ordered state will carry an 
electric charge, giving rise to an anomalous Hall effect, 
different from the "conventional" anomalous Hall effect 
in conventional ferromagnets . It should be noted that 
this phenomenology does not occur in usual condensed 
matter systems because the quantum anomaly is often 
cancelled by the so called fermion doubling effect Q. 

Recently, interaction effects have been investigated in 
the graphene structure, proposing an interesting phase 
diagram based on the quantum Monte Carlo simulation 
[l^ l . Immediately, such a phase diagram was interpreted 
in various analytical scenarios, where the singlet-channel 
interaction was emphasized [l5| . On the other hand, we 
focus on the role of spin-flip scattering, allowed by in- 
teractions in the spin-triplet channel. We introduce an 
anomalous electron self-energy associated with the spin- 
flip scattering, and discuss its possible implication in con- 
nection with the topological nature of the surface state. 



B. Eliashberg approximation 

One can perform an Eliashberg approximation for the 
effective Lagrangian [Eq. (1)] at the quantum critical 
point {m\ = 0), where both fermion and boson self- 
energy corrections are determined fully self-consistently 
in the one-loop level. Considering the cumulant expan- 
sion up to the second order, we can flnd an effective action 
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Sqc = J J ^^^'(■^1(^)[{(<5t - tJ-)Saa' +iVfZ ■ (cT X V)aa'}6{x - x') + ll^a'ix - x')]lpa'ix') 

— Nc,Tjaa'{x — x')Ga'a{x' — x)^ 

+ J (fx J (fx'(4>,{x)[{~d^ - vlV^)S{x - x')5,j + U,j{x - a;')]0j(a;') - n,j{x - x')D,j{x - a;')) 

+ y / d^x I (fx' Dij{x - x')tr{a,G{x ^ x')ajG{x' - x)y (2) 



Y,aa'{x — x') and Ga'a{x' — x) are the electron self- 
energy and Green's function, respectively, where a and a' 
represent spin states, t and I. Ily {x—x') and Dij [x—x') 
are the spin-fluctuation self-energy and Green's function, 
respectively, where i and j express spin directions, x and 
y. All repeated indices are summed. The last term is 
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called the Luttinger-Ward functional 16f], which is the 
key ingredient for an interaction effect obtained in the 
Eliashberg approximation. 

Integrating over electrons and spin fluctuations, we 
obtain an effective free energy as a functional for self- 
energies 



N f (Pk r ~i 

i^w[I]M,n(q,ir!)] = -^^ / --^<w|ln[-G-i(fc,zc^)] + S(fc,zc^)G(fe,zc^)| 

+^-Y^E / (Sfs^E / -^D,,{q,iil)tr^^,(^(T,G{k + q,zuj + iil)a,Gik,ZLj)y (3) 



where 



g ^{k, ioj) — S(fe, 



Gik,iuj) 

D{q,in) = ([n^ +v^q^]I + U{q,in)) ^ (4) respectively. 



are fully renormalized propagators of electrons and spin 
fluctuations, respectively. 



g{k,iuj) 



{iLJ -\- fl)I + VftijkiCTj 



v'jk'^ 



is the bare propagator of electrons, where is the two 
dimensional antisymmetric tensor with i,j — x,y. N„ 
represents the spin degeneracy, in our case = 2. traa' 



and trij mean trace for spin states and spin directions, 



Minimizing the free energy functional with respect 
to both self-energy corrections, we obtain fully self- 
consistent coupled Eliashberg equations for electron and 
boson self-energies 



1 . f (Pk ( 
TLij{q,ift) = Na-^-^'^ J j^—j^traa'[(TiG{k + q,iuj + in)(7jG{k,iuj) 

'E{k,iuj)^g^-J2 J j^D,j{q,in)(T,G{k + q,iuj + in)aj. 



(5) 



As discussed before, the triplet interaction channel gives diagonal term in the electron self-energy matrix. The 
rise to the spin-flip scattering, corresponding to an off- normal Eliashberg self-energy is proposed to depend on 
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frequency only because the momentum dependence is 
regular and the singular behavior can be shown from the 
frequency dependence 0, On the other hand, the 

momentum dependence for the anomalous self-energy has 
not been clarified yet. We propose the following ansatz 
for the electron self-energy 



(6) 



where the momentum dependence for the off-diagonal 
self-energy is given by the same dependence as the topo- 
logical band structure. In other words, the off-diagonal 
self-energy is added to the dispersion relation as the nor- 
mal self-energy in the conventional system, k = k/\k\ is 
the unit vector. 



C. Polarization function 

Inserting the self-energy expression into the Green's 
function, we obtain 



G{k,iuj) 



+ fi — Y^{ibj)]I + [vf\k\ + ^{iijj)]eijki(Tj 
[iij + fi- S(ia;)]2 - [vf\k\ + $(iw)]2 



Based on this expression, we can evaluate the boson self- 
energy, given by the fermion polarization bubble. A de- 
tailed procedure is shown in appendix A. 

It turns out that spin dynamics at the XY ferromag- 
netic quantum critical point is given by the following ef- 
fective Lagrangian 



J 



= (ClO' 



vl\q\'')ct^M,i^)PM-<l^-i^) + 



f]2 



r 



where 



is the projection operator to satisfy PikPkj = Pij- Cj^ct) 
is a constant of the order of 1 and wl(t) is the renormal- 
ized velocity for spin fluctuations, where L{T) represents 
the longitudinal (transverse) mode. An important fea- 
ture for spin dynamics is multi-scale quantum criticality, 
where the longitudinal spin dynamics is given by the dy- 
namical exponent z — 3 while the transverse one is de- 
scribed by z = 2. The multi-scale quantum criticality can 
be interpreted as follows. The longitudinal spin dynam- 
ics is not involved with the spin-flip scattering, described 
by the z = 3 ferromagnetic quantum criticality. On the 
other hand, the transverse dynamics is associated with 
the spin-flip process. In the topological surface state this 
process appears with finite momentum transfer, similar 
to the process with antiferromagnetic fluctuations. As 
a result, the transverse spin dynamics is described by 
z = 2. We note that the self-consistent Eliashberg ap- 
proximation gives essentially the same result as the ran- 
dom phase approximation [6j]. 



where Z?i(T) (9, ii^) is the kernel of the longitudinal 
(transverse) mode. Based on this spin-fluctuation prop- 
agator, we find an electron self-energy in the Eliashberg 
approximation. An important point is the role of the 
multi-scale quantum criticality in the electron self-energy. 
It turns out that the electron self-energy is governed by 
the longitudinal spin dynamics at low energies because 
the z = 3 critical dynamics gives rise to more singular 
corrections than the z = 2 transverse mode. Actually, 
this result was expected in the Eliashberg approximation. 
But, it was claimed that a more complicated structure is 
hidden beyond the Eliashberg approximation ,18i] , which 
will be discussed in the last section. 



D. Electron self-energy 



The spin-fluctuation propagator is given by 



Dij{q,in) 



DL{q,in)P,j 



DT{q,tn){S^j - P^j), 



r n2 



(8) 



Our novel point is that the anomalous self-energy cor- 
rection is exactly the same as the normal one, given by 



5 



where the last approximate equaUty is given by the z = 3 
longitudinal mode DL{q,iil). This equivalence is at the 
heart of the self-consistency in the Nambu-Eliashberg 
framework. As discussed in the introduction, the self- 
energy correction in the fermion bubble diagram changes 
the z — 2 dynamical exponent in the transverse mode. 
Actually, one can see this effect from Eqs. (A9) and 
(A13) in appendix A. If the anomalous self-energy is 
not introduced, the normal self-energy contribution dom- 
inates over the bare frequency, modifying the dynamical 
exponent of the transverse mode from z = 2 to z = 12/5. 
However, introduction of the anomalous self-energy can- 
cels the normal self-energy exactly in the low energy 
limit, recovering the RPA result [6]. 

It is interesting to compare the mechanism of cancella- 
tion with the perturbative analysis in nematic quantum 
criticality, where next leading corrections given by elec- 
tron self-energy and Maki- Thompson vertex corrections 
are introduced for dynamics of critical bosonic modes 
0. Since the fermion self-energy is driven hy z — 3 crit- 
ical fluctuations, the self-energy correction in the trans- 
verse mode, corresponding to the z = 2 dynamics in the 
RPA level, gives rise to an additional singular correction, 
changing the dynamical critical exponent from z = 2 
to z = 12/5. This weird result is cured by the Maki- 
Thompson vertex correction, where the singular correc- 
tion of the self-energy is cancelled by the singular vertex 
correction exactly, which is a fundamental cancellation 
based on the Ward identity. 

The self-consistency of the Nambu-Eliashberg theory is 
far from triviality for multi-scale quantum criticality, and 
this effective field theory is certainly beyond the simple 
Eliashberg description in respect that vertex corrections 
are introduced in some sense. However, it turns out that 
critical exponents for thermodynamics is completely the 
same as the conventional Eliashberg theory without an 
anomalous self-energy correction [19j. 



E. Stability of the Eliashberg framework 

The Eliashberg theory is non-perturbative in respect 
that quantum corrections are introduced fully self- 
consistently in the one-loop level, consistent with one- 
loop renormalization group analysis. Indeed, the scaling 
expression of the free energy was derived explicitly, based 
on the Eliashberg approximation fS]. In this respect the 
Eliashberg theory for quantum criticality implies the the- 
ory with critical exponents given by the Eliashberg ap- 
proximation, satisfying the scaling theory. 

One cautious person may criticize the Eliashberg 
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approximation, neglecting vertex corrections for self- 
energies. Recently, it has been clarified that two di- 
mensional Fermi surface problems are still strongly inter- 
acting even in the large- limit [l9l - l23 |. implying that 
vertex corrections should be incorporated. An impor- 
tant question is whether these vertex corrections give 
rise to novel critical exponents beyond the Eliashberg 
theory. Several perturbative analysis demonstrated that 
although ladder-type vertex corrections do not change 
critical exponents of the Eliashberg theory, Aslamasov- 
Larkin corrections result in modifications for such critical 
exponents ^21, 22]. If this is a general feature beyond this 
level of approximation, various quantum critical phenom- 
ena dl] should be reconsidered because novel anomalous 
exponents in fermion self-energy corrections are expected 
to cause various novel critical exponents for thermody- 
namics, transport, and etc. However, considering our re- 
cent experiences on comparison with various experiments 
in heavy fermion quantum criticality, we are surprised at 
the fact that critical exponents based on the Eliashberg 
theory explain thermodvnamics [2611 . both electrical and 
thermal transport coefficients 12711 . uniform spin suscep- 
tibility 28], and Seebeck effect [23] quite well. 

Recently, one of the authors investigated the role of 
vertex corrections non-perturbatively, summing vertex 
corrections up to an infinite order [l^ [23| • It turns out 
that particular vertex corrections given by ladder dia- 
grams do not change Eliashberg critical exponents at all 
[24| . consistent with the perturbative analysis. This was 
performed in a fully self-consistent way, resorting to the 
Ward identity. In contrast with the previous perturba- 
tive analysis, Aslamasov-Larkin corrections were shown 
not to modify the Eliashberg dynamics [l^ . Resorting to 
the analogy with superconductivity, where the supercon- 
ducting instability described by the Aslamasov-Larkin 
vertex correction is reformulated by the anomalous self- 
energy in the Eliashberg framework of the Nambu spinor 
representation [s^l, we claimed that the off-diagonal 
self-energy associated with the 2kF particle-hole chan- 
nel incorporates the same (Aslamasov-Larkin) class of 
quantum corrections in the Nambu spinor representa- 
tion. We evaluated an anomalous pairing self-energy 
in the Nambu-Eliashberg approximation, which vanishes 
at zero energy but displays the same power law depen- 
dence for frequency as the normal Eliashberg self-energy. 
As a result, even the pairing self-energy correction does 
not modify the Eliashberg dynamics without the Nambu 
spinor representation, resulting in essentially the same 
scaling physics for thermodynamics at quantum critical- 
ity. However, we admit that this issue should be investi- 
gated more sincerely, resorting to a direct summation for 
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such a class of quantum corrections. 



III. ARTIFICIAL ELECTRIC FIELD IN THE 
ENERGY-MOMENTUM SPACE 

An essential question is on the role of the anoma- 
lous self-energy in transport. Considering that the off- 
diagonal self-energy has the same momentum dependence 
as the kinetic-energy term, it is natural to examine its 
role in the Hall conductivity. 

An effective equation of motion for quasi-particle dy- 
namics was proposed as follows [10] 



dR 

~dt 
dk 

H 



(B-Sxv) 



dk 



= -E + B X 



dR 

It' 



(10) 



describing the wave-packet dynamics in the semi-classical 
level. R represents the center coordinate of the wave- 
packet and fe expresses the Bloch wave vector. E and 
B are external electric and magnetic fields while £ and 
B are artificial electric and magnetic fields in the en- 
ergy and momentum space. 3 is usually referred as the 
Berry curvature, v — is the group velocity, where 
Ek is the quasi-particle energy dispersion. A key ingre- 
dient is that interactions give rise to an artificial elec- 
tric field in the energy-momentum space, modifying the 
quasi-particle dynamics according to the Maxwell equa- 
tion. As a result, the Hall conductivity should be cor- 
rected in the following way 10] 



d^k 



f{Ek){B,~[vy£,,-v^£y\^ill) 



where the Berry curvature is shifted hy £ x v due to 
interactions. f{Ek) is the Fermi-Dirac distribution func- 
tion. In this section we calculate this quantity with an 
introduction of the anomalous self-energy correction. 

We consider the following effective Hamiltonian 
Hf.ff = N{k, U!)-(T, which incorporates the effect of inter- 
actions through the self-energy. The pseudospin vector 
is given by 



N{k,uj) = d{k) + Ad{uj), 

d{k) = -Vfky&x + Vfkx&y + dz{k)&z, 

Ad{uj) = -^{uj){ky/k)d-x + <f>(ci;)(fc^/fc)o-y 



The Berry curvature and artificial electric field are ex- 
pressed in terms of the pseudospin vector 

B,^^dk^Nxdk^N-N, 
£y^^d^Nxdk^N-N, 



where the d{k) vector is nothing but the bare dispersion 
and Ad{Ld) is associated with the anomalous self-energy 
correction. dz{k) is introduced for time reversal symme- 
try breaking as the z-directional magnetic field, set to be 
dz{k) — > m. Adz{io) is also introduced for consistency. 

k — ^/kl + ky is the amplitude of the momentum. 



(13) 



where TV = AT/] AT] is the unit vector. 

One can understand this expression based on another 
equivalent representation, where the field strength is ex- 
pressed by the gauge field. Resorting to the CP^ repre- 
sentation [sil 



N-(T^ ]iV]J7(T3C/t, 



(14) 



where U is an SU(2) matrix field to describe a direction 
of the pseudospin vector, we can define an SU(2) gauge 
connection 



(15) 



This naturally leads to both Berry curvature and artifi- 
cial electric field 



£x = dkyA^ - d^Aky, 



(16) 



where the U(l) projected component of the SU(2) gauge 
field is given by A^, = A^^'^ . 

Inserting the pseudospin vector Eq. (12) into Eq. (13), 
we obtain both Berry curvature and artificial electric field 



1 [f/-f $(w)/fc]2[m + Ad,(u;)] 



(^Vfk + $(w)]2 + [m + Adz{uj)]- 



3/2 ' 



1 



where T is given by 

[9^$(c^)](l/fc)[i;/ + $(c^)/fc][m + Adz{uj)\ 



{vjk + $(w)]2 + [m + A4(a;)]2 
[«/ + $(w)/fc]2[a„A4(w)] 



3/2 



Adz{yi)crz, 
(12) 



{^Vfk + $(cj)]2 + [m + A4(cj)]2) 
Considering the group velocity 



3/2 ■ 



Vy — fcyV, — kxV, 



with 



[vfk + <i>{uj)]vf 



k ^[vfk + $(w)]2 + [m + A4(^)]2 ' 
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we find that the contribution of the electric field £ to the 
shift of B in Eq. (11) vanishes identically 

On the other hand, the Berry curvature should be mod- 



ified in the presence of the self-energy correction. As a 
result, we obtain the following expression for the Hall 
conductivity 



-if! [°° [f/fc + $(^)]^[m + A4(^)] _le^ 



where 

F[kH oo] = 

F[kL ^ 0] = - 



[ra + J\dz{(jj)\^{(jj) [m + Mz{uj)]^'^{io) 



(F[kH ^ oo] - F[kL ^ 0] 



(17) 



■In 



{[m + A4(^)]2 + ci>2(^)) ^^2^^^ ^ ^ A4(c.)]2)'^' + yi^RTFTAdlMP 

[m + l\d,Sioyf + [m + A4(w)]$2(a;) 



($2(^) + [m + A4(w)]- 
[m + Ad,(a;)]$2(t^) 



3/2 



(ci,2(^) + [^ + A4(^)]2)'/' [m -t- A4(^)]2 + $2(0.) + V<i>2(c.) + [m + Ad,(a.)]2^ci>2(,.) + + A4(a.)]^ 



It is straightforward to see that the noninteracting case 
in $(w) — recovers the well known result, that is, the 
half Hall conductivity where F[kH — > oo] = and 
F[kL 0] = -sgn{m). 

An interesting point in this expression is that the Hall 
conductivity diverges in the fc^ limit, seen from the 
log term. Although we do not understand the reason for 
this divergence, it implies that we should incorporate ver- 
tex corrections even for the intrinsic topological effect in 
the Hall conductivity when interactions are introduced. 
If we take into account vertex corrections, this log di- 
vergence may be re-summed to result in the common 
power-law behavior with an exponent, which vanishes in 
the fci limit. 

We suspect that the underlying mechanism for this 
divergence may be charged vortices. As discussed be- 
fore, the parity anomaly assigns an electric charge to 
a vortex excitation 131. From the ordered side in the 



duality picture, such vortex excitations are well defined 
particles, thus taken into account with electrons on an 
equal footing. An interplay between electrons and vor- 
tices would be described by the mutual Chern-Simons 
theory 3] , where their mutual statistics is guaranteed by 
the mutual Chern-Simons term. It is not clear at all how 
such an interplay affects transport, in particular, the Hall 
conductivity. It will be an interesting future direction to 
investigate this problem in the duality picture because 
topological excitations carry nontrivial fermion quantum 
numbers in this situation. 



IV. DISCUSSION AND SUMMARY 

In this study we constructed an Eliashberg framework 
for the Hertz-Moriya-Millis theory 01, describing ferro- 
magnetic quantum criticality in the surface state of three 
dimensional topological insulators. An idea is to intro- 
duce an anomalous self-energy correction beyond the pre- 
vious study 0], resulting from spin-flip scattering. We 
could show that our ansatz for the anomalous self-energy 
allows a set of fully self-consistent solutions for electrons 
and spin fluctuations. 

The spin-fluctuation dynamics was shown to coincide 
with the solution based on the random phase approxima- 
tion, where both 2 = 3 and z — 2 quantum critical dy- 
namics arise f6^ . The 2 = 3 quantum criticality describes 
the longitudinal spin dynamics, where only forward scat- 
tering contributions result in the z = 3 Landau damp- 
ing, while the z = 2 quantum criticality for transverse 
spin fluctuations originates from the spin-flip scattering, 
involved with the finite-momentum transfer due to the 
spin-orbital quenching. 

This multi-scale quantum criticality was shown to al- 
low the self-consistent Nambu-Eliashberg solution for the 
electron dynamics. The normal self-energy turns out to 
follow the 2 = 3 quantum criticality because such critical 
fluctuations give rise to more singular contributions than 
the 2 = 2 transverse mode. We found that the anomalous 
self-energy shows exactly the same frequency dependence 
as the normal self-energy although the underlying mech- 
anism is not completely clear. The equivalence between 
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the normal and anomalous self-energies with a different 
sign is at the heart of the self-consistency in the Nambu- 
Eliashberg theory for multi-scale quantum criticality. In 
other words, the simple Eliashberg theory without the off 
diagonal self-energy cannot be self-consistent for multi- 
scale quantum criticality. We argued that the mechanism 
of cancellation between the normal and anomalous self- 
energies in the fermion bubble diagram may be rooted in 
the Ward identity, comparing our result with the pertur- 
bative analysis for nematic quantum criticality, where the 
singular self-energy correction is cancelled by the vertex 
correction based on the Ward identity 'oj- 

We demonstrated that such an off-diagonal self-energy 
gives rise to an artificial electric field in the energy- 
momentum space beyond the static Berry curvature from 
the topological band structure. In particular, we calcu- 
lated the anomalous Hall conductivity in terms of both 
Berry curvature and artificial electric field. Although 
the contribution from the artificial electric field vanishes 
identically, we could find that the Berry curvature is 
modified due to the anomalous self-energy correction, 
resulting in the log-divergence for the Hall conductiv- 
ity. We argued that the log-divergence may disappear 
if vertex corrections are re-summed up to an infinite or- 
der, which is expected to turn the log-divergence into the 
power-law dependence, vanishing algebraically. We sug- 
gested an idea that this phenomenology may be related 
with charged vortices in the duality picture, where the 
electron quantum number of the vortex results from the 
quantum anomaly, an essential feature of the topological 
surface state. 

Before closing, we discuss the stability of the Eliash- 
berg approximation in another aspect, considering an 
interesting recent study for nematic quantum criticality 
[l8| . The nematic quantum critical point is well known 
to show multi-scale quantum criticality, where critical dy- 
namics of nematic fiuctuations is governed by both z = 3 
and z = 2. An important notice is that although the 
self-energy correction results from scattering with z = 3 
critical fluctuations dominantly, scattering with z — 2 
critical modes gives rise to the logarithmic singularity 
for the quasi-particle residue pisj , which means that ver- 



tex corrections should be introduced up to an infinite 
order. Based on a complicated renormalization group ar- 
gument, the authors proposed an interesting expression 
for the electron Green's function 

G(k,iuj) cx 7: ° — . 

Do has a constant of an energy unit proportional to the 
Fermi energy, and T,{iuj) oc Iwp/"^ is the z — 3 self-energy. 
Efc is the band dispersion. It should be noted that the 
anomalous exponent rj results from scattering with z ~ 2 
critical fiuctuations. 

In the present system a significant quantity is the off 
diagonal self-energy. Although such an anomalous self- 
energy will not be allowed in nematic quantum criticality 
due to the absence of the spin-orbit coupling, it may play 
an important role for the exponent rj in ferromagnetic 
quantum criticality of the topological surface. A more 
complete framework would be to introduce vertex correc- 
tions fully self-consistently with self-energy corrections in 
the Nambu-spinor representation, where the off-diagonal 
self-energy contribution is also included. 

In summary, quantum criticality with the topological 
band structure is expected to open a novel direction of 
research in respect that the interplay between interac- 
tion and topology may cause unexpected non-Fermi liq- 
uid physics. We found an anomalous behavior in the 
intrinsic topological contribution of the Hall coefficient, 
the source of which is the energy-momentum-space field 
strength beyond the static Berry curvature as a result of 
the interplay between quantum criticality and topological 
structure. 
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Appendix A: Polarization function 



In appendix A we derive the self-energy correction for spin fiuctuations, where the longitudinal spin dynamics is 
given by z = 3 while the transverse one is described by z = 2. Inserting the electron Green's function with the ansatz 
for the electron self-energy matrix into the Eliashberg equation for the boson self-energy, we obtain 



(T 



.9^ 1 [ (fk^ f ^ [iij + iQ + ^ - T,{iu} + + [vf\k + q\ + ^{iuj + iri)]enm{kn + qn)o'7i 

in + fi~ E(ia; + ifl)]^ - [vf\k + q\ + ^{iuj + iil)]'^ 

' Huj + n~ I](iw)]2 - [vf\k\ + $(iw)]2 }' ^^^^ 



[iuj + + fi~ 'E{iuj + «f7)]2 - [vf\k + q\ + ^{iuj + zfi)]^ 



where the following approximation for the momentum 

k 4-n - ^n + Qn kn Qn , 

\k + q\ Ikp] \kF\ 

is utilized. 

Resorting to identities for Pauli matrix products 
<Ti<Tj = SijI + ieijiCTi, 

CiCTmCfjCrm' = Sim(^Sjm' I + iijm'l' Cfl'^ + iiimjC^m' — {SijSms — SisSmj)(^Ssm,' I + i^sm'l'Cfl'^ (A3) 

and performing the decomposition for the denominator, we can rewrite Eq. (Al) as follows 

)- ^9 J (2;r)2 2[iuj + in + - S(iw + ifl)] 2[iuj + ^ - 

+ zf^,w/|fc + q\ + ^{ito + in)]g[iLj,Vf\k\ + $(icj)] + g[iuj + iQ, -Vf\k + q\ - ^{iuj + in)]g[iuj,Vf\k\ + 
+g[iij + iQ.,Vf\k + q\+ ^{ico + ifl)]g[iLo, -Vf\k\ - $(iw)] + g[iu + in., -Vf\k + q\- ^{ico + in)]g[iuj, -Vf\k\ - 

+ + - E{iuj + iQ.)][iuj + M - - [vf\k + q\ + ^{iuj + iQ.)][vf\k\ + ^{iLj)]{k + q) ■ k^dij 

+ [vf\k + q\ + + iO)]['!;/|fe| + $(ia;)][ej„ej„'(fe„ + q„)fe„' + ejnein'{kn + 9n)fen']}, (A4) 

where is an antisymmetric tensor with i,j = x, y, and 

I 

g[iu), X] 



icj + fj, — Y,{ico) — X 
We rearrange the above expression as 

d^k „ 



u,jiq,m) = nJ-^Y1 J ■ 



(27r)2"'^ 

^^^g[iuj + iCt,Vf\k + q\ + ^{ico + iQ.)]g[ico,Vf\k\ + ^{icj)] + g[ico + iO., -Vf\k + q\ - ^{iuj + in)]g[iuj,Vf\k\ +^{iui)] 

+g[iij + in,Vf\k + q\ + ^{iij + iQ)]g[iij, -Vf\k\ - + g[ico + i^, -Vf\k + q\ - ^{ioj + i0,)]g[iij, -Vf\k\ - 

1 . f cPk ( - - ^ ^ - ,1 

+^<^^^'^ j I ' ^^^^ ^ l^tnejn'iK + qn)kn' + ijn<^tn'{kn + qn)fcn']| 

^g[ii^ + iQ.,Vf\k + q\ + ^{iuj + iQ)]g[iuj,Vf\k\ + - g[iuj + iQ-,-Vf\k + q\ - ^{ioj + iVi)\g[ii^,Vf\k\ +^{iuj)] 

-g[ii^ + iQ-,Vf\k + q\ + ^{ii^ + iVl)\g[iu,-Vf\k\ - +g[iuj + iQ-,-Vf\k + q\ - ^{iuj + iVi)\g[iu,-Vf\k\ - 

(A5) 

where the frequency term in the numerator is reorganized. 

It is convenient for the momentum integration to linearize the band dispersion near the chemical potential 

— Vf\k + q| ~ —Vfkf ■ (fe — kf) — Vfkf ■ q, (A6) 
where the Fermi momentum kf = \kf\kf is defined from /z — Vf\kf \ = 0. 
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Inserting the linearized band into Eq. (A5) and performing the decomposition for the denominator, we obtain 
,g2 1 ^ f (Pk ^ 1 



(27r)2 ifi _ + ifl) - S(iw)] - [^{icj + iQ) - ^{icj)] - Vfkf ■ q 

I . . 

'^iuj — — <E>(icj) — Vfkf ■ {k — kf) ioj + ifl — S(za; + iQ) — ^{ioj + ifl) — Vfkf ■ {k — kf) — Vfkf ■ q J 

1 . f cPk r ^ -~ ^ ~ 



iQ - [E{ioj + ifl) - E(iw)] - [^{ioj + ifl) - ^{ioj)] - Vfkf ■ q 

I ^ ^ ^ U7) 

'-ioj — T,{iu)) — ^{iuj) — Vfkf ■ {k — kf) ioj + iCl — T,{ico + ifi) — ^{iui + iQ) — Vfkf ■ {k — kf) — Vfkf ■ q ' 

where only dominant contributions are selected, which means that all terms with the chemical potential in the 
denominator are safely neglected. 

Before evaluating integrals, we summarize each component for the spin-fluctuation self-energy as 



^ iui — 



(Pk 1 



r)2 in - [S(«w + in) - - [<^{iuj + iVi) - - Vfkf ■ q 



1 1 



iu3 — Yj{iuj) — — Vfkf ■ (fe — kf) iu) + — S(ia; + ifi) — $(iw + iO) — Vfkf ■ (fe — kf) — Vfkf ■ q ■ 
g^lsr f -{k + q) ■k + 2{ky + qy)ky 



4 ^ ^ J (27r)2 iQ _ [S(ia; + iQ) - T,{iuj)] - + iO) - $(ia;)] - Vfkf ■ q 

I ! I ^ I ^ _\ 

iw — E(iu;) — $(ia;) - vfkf ■ {k - kf) iw + ifi — Y,{iuj + ifi) — $(iw + ifi) - Vfkf ■ {k — kf) — Vfkf ■ q J 



XT / -r^x AT 9 1 \ - [ d^k 
n,,(,,^fl) = Ar.^-^y^ 



. . ,^ ifi - [I](ia; + if^) - E(iw)] - [$(iw + if^) - i>(iw)] - -u/fe/ • q 

I ^ ^ —I 

iw — S(ia;) — $(iw) — Vfkf ■ (fe — fc/) iw + ifi — T,{ioj + ifi) — $(iw + ifi) — u/fe/ • (fe — fe/) — Vfkf ■ q ' 

+N ^1 V f — -(fc + g)-fc + 2(fc^+g:,)fc^ 

" ^ I^^J (2'r)2 in - [E{iui + in) - E(iw)] - [$(iw + if^) - $(iw)] -Vfkf-q 

I ^ ^ ^ ^ _\ 

iw — S(ia;) — ^{icu) — Vfkf ■ (fc — kf) iw + ifi — S(ia; + ifi) — $(ia; + ifi) — Vfkf ■ (fc — kf) — Vfkf ■ q ' 

n (ozf!) = n (a in) = -N ^-S^ f ^ 2fc + (gj, + g,fc.) 

xyyi, ) yxKH, ) ^ 4. ji J {2'k)'^ iQ - [Y.{iLO + iQ) -^{iuj)] - [^{iuj + in) - ^{iw)] - Vfkf ■ q 

I ^ ^ ^ ^ ^ i(.A8) 

iw — E(iu;) — <^(iuj) — Vfkf ■ (fe — fe/) iw + ifi — E(ia; + in) — ^(iuj + ifi) — Vfkf ■ (fe — fe/) — I'/fe/ • q 

First, we consider the diagonal component of the boson self-energy. It is convenient for the momentum integration 

to take the angular coordinate 

^xx{q,i^) 

-"^ 1 '•^^ .„ r . 2sin2 6' 



_ ^ 1 n 2sin^g 

~ " 'ivf{2n)'^ /3 ^ Jo ^ in - [E{ioj + in) -T,{iw)]-[^{icv + in) -^{icv)]-Vfq:c cose -VfQy sine 

I io 



1 1 



E(ia;) — <l>(icj) — e iw + ifi — E(ia; + ifi) — ^{iuj + in) — e — Vfq^ cos 9 — vfQy sin0 J 

= N ^^g' r — r de 2sm^9[sgn{Lj)-sgniu: + n)] 

''4w/(27r)2 27r io if^ - [E(icj + if7) - E(iw)] - [$(iw + ifi) - $(iw)] - u/fe cos6' - u/g^ sin6'' ^ ' 

where we considered the zero temperature limit in the last equality. 
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The angular integration is given by 







zil' — Vfqx cos 9 — VfQy sin( 



where 



n' = n + i[j:{ioj + in) - I](ia;)] + i[<S>{iuj + in) - ^{icj)] 



is an effective frequency with both normal and anomalous self-energies. 

Using this angular integration, we find that both z = 3 and 2 = 2 critical dynamics appear in the following way 

2 roc 

Il.j,x{q,in) oi N^- — / —[sgn{uj) - sgn{uj + n)] 

{[ 1 ^ n ^ n^ 1 ( ^ ^ ^ yh 

IV /ri2 I „,2„2 l'?0'^ W?a^ /r>2 I „.2„2/ 0^ \W?o2 V?o2 /r)2 I „,2„2 / a-^ J 



+ W^gZ 



^N.^o^^m^^y4-^^^m^^)i (All) 

A key point in our derivation is fi' w 17 in the low energy limit. In other words, the normal self-energy is cancelled by 
the anomalous self-energy exactly in the low energy limit, explicitly shown in appendix B. If we do not introduce the 
anomalous self-energy correction, the normal self-energy dominates over the bare frequency, changing the dynamical 
exponent from 2 = 2 to z = 12/5 for the transverse mode, consistent with the perturbative analysis in nematic 
quantum criticality @. 

The other diagonal boson self-energy can be obtained in the similar way 

8TTVf\Vfq vjq^ v^q^ J q^ SnvfKvjq^ vjq^ J q^ 

Rewriting the off-diagonal boson self-energy in the angular coordinate 

2 cos 6 sin 9 



in ~ [I](ia; + in) — I](za;)] — [<I>(iw + in) — — Vfqx cos9 — Vfqy sin9 

/ I . I _\ (A13) 

liuj — I]{iuj) — $(ia;) — e iuj + in — ^{iuj + in) — $(ia; + in) — e — v/qx cos 9 — Vfqy sin 9 i ' 

we reach the final expression 

STTVf \Vfq vjq'' v^q'^ / (7^ 

It is straightforward to rewrite Eqs. (All), (A12), and (A14) with an introduction of the longitudinal projection 
operator Pij. Then, we find the 2 = 3 critical dynamics for the longitudinal mode and the 2 = 2 critical dynamics 
for the transverse mode, given by Eq. (7). 



Appendix B: Electron self-energy 

In appendix B we derive an electron self-energy. In particular, we show that the ansatz of Eq. (6) allows the fully 
self-consistent Eliashberg solution, where the off-diagonal electron self-energy turns out to be the same as the normal 
self-energy. Inserting the electron Green's function into the Eliashberg equation for the electron self-energy, we obtain 

. , 2j_V^ f d'^Q n / ■o\ + in + ^ - J:{iuj + in)]I + [vf\k + q\ + ^(iui + in)]enm{kn + qn)crrn 

^ "^''^j^J {27r)^ [tLo + tn + fi-i:{tio + in)]^-[vf\k + q\ + ^ioj + tn)]^ 

(Bl) 
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Considering Pauli matrix identities Eq. (A3) and taking the decomposition for the denominator, we can separate 
the electron self-energy matrix into two parts 

S(fe,ia;) = Si(fe,ia;) + S2(fe,iw), (B2) 

where each part is given by 

Si(fc,zu;) = ^i^ I -^(DI^{q,^n)+DT{q,^n)) 

f 1 1 ^j. 

Uw + ifj + /i - Y,{iuj + ?:0) - vj\k + q\ - + iO) iuj + iO. + ji - + iVL) + Vf\k + q\ + + iil) / 

-\jY.} T^,^ ^TiFTi >-{Da,,i^) + Dr(,,in)) 

I \ I I 

liuj + + jji- S(iw + iO) -Vf\k + q\- <^{iuj + iO) ioj + iQ. + ji- S(ia; + iQ) +Vf\k + q\+ ^{iu) + ifi) J ' 



^2 



\k + q\ 

r \ \ 

Uw + + /u - ll{iuj + iO) —Vf\k + q\- $(ia; + ifi) + + /x - S(it<; + iO) + i;/|fe + q| + + iO) V 

respectively. We will see that S2(fe,iw) is irrelevant at low energies, compared with Si(fc,«a;). 

First, we consider Performing the linearization near the chemical potential [Eq. (A6)] and taking only 

dominant contributions near the Fermi surface, we obtain the following expression 



5^1y- f ^ DL{q.i^)+DT{q.i^) ^ 

2/3^7 (27r)2 iuj + in-T.{iu + in) - <P{iuj + iQ) - Vfkf ■ {k - kf)-Vfkf q 



-2 1^/- (fq /i- _ ,^ _ ^ DL{q,i^) + DT{q,i^) 

+ iO — T,{iu + ifi) — + ifi) — Vfkf ■ [k — kf) — Vfkf ■ q 



2/3^i (27r)2VfcF ^ A;f + ifi - E(ia; + ill) - + ifi) - t;/fe/ • (fe - fe/) - t;/fe/ • q' 

Introducing the polar coordinate 

q ~ q (cos 6x + sin 6y) , 

fe/ = coS(^!E + sin^y, (B5) 
we can rewrite Si(fe,ia;) as I]i(fe,ia;) w E(ia;)J + ^{iuj)eijkiaj + ASi(ia;,feF), where 

ASi(ia;,feir) = V / dgg / (/6»( cos(6i + ?!.)o-j/ - sin(6l + (A)<Tx ) 
£>i(g,iO) + DT(9,iO) 



+ iO — S(ia; + ifi) — ^{icj + iQ) —Vfq cos 6 



(B6) 



can be regarded as a correction for the Eliashberg self- energies, E(iw) and $(iw) given by Eq. (9). 



It is straightforward to perform the angular integral, reaching the following expression 
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9^ 1 f°° q ( \ ^ » 

ASi(iw,fcF) = ;^ — oVl / dqq—\DL{q,i^) + DT{q,i^)]{kxCTy -ky(Tx) 
IVfkp P ^ Jo I^F ^ ' 



r dqq^(^DL{q,in) + DT{q,in)) 



[w + O + iT,{iuj + i^l) + i^{ioj + 2fi)]sgn(w + O) 



2vfk^ 



2vfki 



a/[w + + zS(iw + in) + + iQ)]'^ + {vfqf 

iijj + in — Y^{iijO + in) — <^{i(jj + in) — Vfq\f_ j- 

(B7) 



Comparing this expression with the Eliashberg self-energies, we find that ASi(iw,fcF) is irrelevant at low energies 
due to the fact that the scaUng dimension of the argument itself in the integral expression is higher than that of the 
Eliashberg solution in addition to the additional momentum factor q/kp- As a result, we obtain 



Si(iaj,feir) = T,{iuj)I + '^{iijj){k,c(Ty — kyCTx), 



'£{iu>) = -^ioj) 



^ in 



dqq(^DL{q,in) + DT{q,in) 



sgn(a; + n) 



^[w + n + iT,{iu! + in) + i^{iu! + in)]'^ + vjq^ 

(B8) 



Next, we show that S2(iw, kp) is irrelevant at low energies. Performing the linearization near the chemical potential 
[Eq. (A6)] and taking only dominant contributions near the Fermi surface, we obtain the following expression 



^2{iuJ,kF)^g^^Y.j 7S2 { - (^^ (9' 4 + (g, if^) 4) fej/^x + (Dl {q, in) ^ + DT{q,in)^) kx^y } 

iQ, ^ ^ ^ ^ ^ ^ 



^a; + iO — T,{iuj + iQ) — ^(iuj + iO) — Vfkf • q 



^a; + iO — S(icj + iO) — ^{iu + iQ.) — Vfkf • q 
1 f (fi 

J J^{[DUq,tn) - DT{q,in)]^) {kxCTx - ky^y) - 
/ ^{[DL{q,in)-Driq,in)f''''^^'^^ 



ioj + in — J^{ioj + in) — ^{ioj + in) — Vfkf ■ q 
1 



q^ J\kF kp iui + in-i:{iio + in) - ^{iui + in) - Vfkf ■ q 

(B9) 
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Introducing the polar coordinate of Eq. (B5) , we rewrite the above expression as fohows 

1 f°° f^^ r / \ 
S2(ia;,feF) = -^-y^ / dqq ddl -(DL(q,i^l) coa'^ 9 + Driq^i^l) sm"^ 9] sin (pa 
^^^mJ^ Jo ^ ^ ) 

+ [Dl {q, in) sin^ 9 + Driq, ^^^) cos^ 9) cos4><Ty | . ^ ^ 1,. ^ 7 

V / ) lu + III ~L,(iuj + III) — fvyiuj + lit) — Vf cos{ 

„2 1 roo n2-K 

+^-^y dggy d(?|-(^i?L(g,«f^)cos2f? + ^T(q,«f^)sin2(?j A. singer. 



i?i(q, zr!) sin2 + DT[q. t^) cos^ 0) cosg^ J - ^ -n vr ^ 'OA d>r ^ 1 

^-^y dqgy d6'[L»L(q,ir2) - Z?T(g, if^)](^cos0(T^ - sin^o-j^j 



zo; + — S(ia; + if2) — $(za; + iJl) — f / cos( 

r,2 1 roo n2-K 



J-^^'^ j dqq j d9[DL{q,ii^) - DT{q,in)]^(^-^cos9(Tx - ^sm9cry 



iuj + iQ — + ifi) — ^{iuj + ifl) — Vf cos( 



(BIO) 



Performing the angular integration, one can see irrelevance of this contribution, basically resulting from higher order 
momentum integrals. We prove that our ansatz is fully self- consistent in the Nambu-Eliashberg approximation. 
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